Electronic Green's functions in a T-shaped multi-quantum dot system 
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We developed a set of equations to calculate the electronic Green's functions in a T-shaped 
mult i- quantum dot system using the equation of motion method. We model the system using a 
generalized Anderson Hamiltonian which accounts for finite intradot on-site Coulomb interaction 
in all component dots as well as for the interdot electron tunneling between adjacent quantum 
dots. Our results are obtained within and beyond the Hartree-Fock approximation and provide a 
path to evaluate all the electronic correlations in the multi-quantum dot system in the Coulomb 
blockade regime. Both approximations provide information on the physical effects related to the 
finite intradot on-site Coulomb interaction. As a particular example for our generalized results, we 
considered the simplest T-shaped system consisting of two dots and proved that our approximation 
introduces important corrections in the detector and side dots Green's functions, and implicitly in 
the evaluation of the system's transport properties. The multi-quantum dot T-shaped setup may 
be of interest for the practical realization of qubit states in quantum dots systems. 

PACS numbers: 73.63.Kv,72.15.Qm,72.10.-d 



I. INTRODUCTION 

In the last decade, quantum dot (QD) systems pro- 
vided the perfect environment for the study of com- 
plex many-body effects such as the Kondo effect and the 
Coulomb blockade For example, electronic transport 
throughout single QD systems allow the controlled real- 
ization of the Kondo regime of the Anderson single impu- 
rity problem^. On the other hand, electronic transport 
in more complex systems, such as double QD arranged 
in a series, parallel, or T-shaped configuration Q, can 
be explained based on the interplay between Kondo res- 
onances and Fano interference effects These sophis- 
ticated physical phenomena can be precisely controlled 
in mesoscopic devices using external parameters such as 
bias voltages or external electromagnetic fields. 

Theoretically, a simple model to account for the sys- 
tem's QD's is the Anderson model @, each dot being 
represented by a localized level similar to the impurities 
in the original problem. Additional terms are introduced 
to describe the interactions between the system's QD's 
or the interactions between the QD's and the external 
leads required by the transport studies. Of main impor- 
tance is the on-site Coulomb interaction term, which reg- 
ulate the occupation number in each QD, leading to the 
direct expression of the Coulomb blockade in transport 
phenomena^]. Consider the case of a T-shaped double 
QD system, with one dot (detector dot) directly con- 
nected to the external leads and the second dot (side 
dot) coupled to the first one but not to the external leads. 
Various approximations were applied to estimate the role 
of the on-site Coulomb interaction in this case. Wu et 
al. 0] considered an infinite on-site Coulomb interaction 
in the detector dot so the double occupancy was forbid- 
den in this dot and neglected it in the side dot. On the 
other hand, Guclu et al. I8| assumed an infinite on-site 



Coulomb interaction in the side QD, neglecting it in the 
detector QD. Tanaka and Kawakami [3| considered the 
on-site Coulomb interaction to be infinite in both compo- 
nent QD's. Although these configurations are different, 
both the Kondo effect and the Fano interference effect 
play an important role in the system's transport proper- 
ties [1, 0, Hf ■ For each of these studies the main ingredient 
is the evaluation of the system's QD's Green's functions, 
as the transport properties depend directly on the detec- 
tor's dot density of states. 

QD systems were also studied in the nonequilibrium 
regime [3-11 1|. Nonequilibrium transport studies require 
special many body techniques based on the Keldysh for- 
malism such as the perturbation theory and the pertur- 
bative renormalization group fl2rll5j . slave boson[l6lfl7|. 
or the equation of motion (EOM)[l8[. Additionally, 
several studies used numerical methods such as the 
Numerical Renormalization Groupfl9j or the Quantum 
Monte Carlo [20j. Although each of these approaches pro- 
vided reasonable answers for the transport properties in 
nonequilibrium mesoscopic systems, all approximations 
have their limitations. Again, the evaluation of the elec- 
tronic Green's functions in each of the component QD is 
the key ingredient of the calculation. 

Here we propose an investigation of the T-shaped 
multi-QD's system (See Fig. [J) based on the EOM 
method. The system consists on N quantum dots; the 
detector dot (1) is connected to the external leads R 
and L, while the rest of the dots are connected in a 
chain arrangement to the detector dot. Our calcula- 
tion is performed for the most general case when the 
on-site Coulomb interaction has a finite value in each of 
the system's QD's. We consider two different approxi- 
mations for the evaluation of the system's Green's func- 
tions. First, we follow the method introduced by Hew- 
son [2l| for the study of the Anderson's single impurity 
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FIG. 1: Schematic representation of the T-shape quantum dot 
system. The detector dot (characteristic energy ei) is coupled 
both to the side quantum dots and the external electrodes L 
and R. 

model and we obtained results similar to the Hartree- 
Fock approximation. Second, we consider an approxi- 
mation beyond the Hartree-Fock approximation which 
allows us to include higher order effects in the calcula- 
tion of the system's Green's functions. We derive general 
recurrence relations which require self-consistent calcula- 
tions for the evaluation of the electronic Green's function. 
As an example we will study the double dot system, and 
prove that the general relations obtained for the gen- 
eral T-shaped multi-QD's system are leading to impor- 
tant corrections when one analyzes the system's physical 
properties. The paper is organized as follows: In Sec- 
tion II we present the general Hamiltonian of the system 
and the general equations for the system's Green's func- 
tions using the EOM method. In Section III we analyze 
these equations and present their solutions within and 
beyond the Hartre-Fock approximation. Finally, Section 
IV presents our conclusions. Details of our calculations 
are provided in the Appendix. 

II. MODEL 

The multi-quantum dot T-shaped system is modeled 
by a generalized Anderson model: 

N 

h = £k4 CT;a c kCT;Q + £ iJ2 a ^ a ^ 

k,cr;a i— 1 o 

N 

+ Y Uin ia ni^ a + Y ^kl;a (cia;a a l<? + a L c k<r;a) 
i— 1 k(T;a 
N-l 

+ Y Y (4x a *+^ + a !+i<r a ^) • (i) 

i=l <r 

The first term in the Hamiltonian describes the free elec- 
trons in the leads, cj^^ and Ck<T ;Q being the fermionic cre- 
ation and annihilation operators for electrons with mo- 
mentum k and spin a in the lead a (a = left (L), right 



(R)). The following two terms describe the mesoscopic 
part of the Hamiltonian and correspond to the electrons 
localized in the detector (i — 1) and side (i = 2, . . . , N) 
QD's of the system; here e% is the energy of the localized 
level and a\ a and ai a are fermionic creation and anni- 
hilation operators for localized electrons with spin a in 
the i'th QD . Additionally, electrons in each component 
dot are subject to finite on-site Coulomb interaction de- 
scribed by the interaction terms Ui\ in this interaction 
term, n ia = a\ a ai<j, is the number of particle operator 
corresponding to the electronic level £j and electron spin 
a. The last two terms in the Hamiltonian describe in- 
teractions between the system's electrons. The coupling 
constant tj.j+i characterizes the electron tunneling be- 
tween two adjacent dots i and i + 1. Vki ;a characterizes 
the interaction between the free electrons in the lead a 
and the localized electrons in the detector dot (i = 1). 
For simplicity we will consider the case Vki ; L = Vki ; i? in 
which the detector couples to the leads only in the sym- 
metric combination c-^ a = (ckcr ; L + Ckcr;i?)/"\/2 and the dot 
connects effectively to a single lead, with Vki = V^Vki-L- 
The system's physical properties can be investigated 
using the Green's function formalism. In particular, to 
investigate the transport properties of the system, one 
will be interested in calculating the detector's dot (i = 1) 
localized electrons Green's function. One way to extract 
the characteristic Green's function for localized electrons 
in the system's QD's is to use the EOM method. It is well 
known that in the case of a general Anderson impurity 
model, the EOM method leads to an infinite hierarchy of 
higher-order Green's functions, so in order to obtain the 
QD's electronic Green's functions one needs to introduce 
a reliable approximation to truncate this hierarchy. The 
difficulty is mainly introduced by the interaction terms in 
the system's Hamiltonian. When the on-site Coulomb in- 
teraction term is absent, an exact solution of the problem 
is possible as it is well known that in this case the set of 
equations obtained from the EOM method are closed. In 
the case of two fermionic operators A and B the Fourier 
transform of the Green's function with respect to the 
time, Gab{u) = ((A;B)), is given by the general equa- 
tion 

u((A;B)) = ({A,B}) + (([A,ir\;B)) , (2) 

where (A) represents the mean value of the operator A, 
{A, B} the anti-commutator of the operators A and B, 
and [A, B] their commutator. The last term on the right 
hand side of the equation is responsible for the generation 
of the infinite chain of higher order Green's functions. 

We start from the evaluation of the electronic Green's 
functions corresponding to the conduction electrons, 
^kk'C 1 ^) — << c k<r; cj^CT >>, and for electrons localized in 
quantum dot i, Gf^ui) =<< ai a ]a\ a >>. Consider first 
the free electrons Green's function. Using the general 
equation, after some simple calculations we find 

(w - e k ) Gkk' M - V kl G^, (w) - <W , (3) 
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where <5kk' is the Kronecker symbol. On the other hand, 
the localized electrons Green's function will be the solu- 



tion of the following equation: 
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(w - et) G1M - E/inj(w) - Y, ViiG?£i(w) - 6(iV - 1 - i)t M+1 Gf +l! » - 0(i - 2)t i _ 1 , i Gf_ 1 » = 1 , (4) 



r 



where 0(x) is the standard theta function with ©(as) = 1 
for x > and ©(a;) = in rest. One can see from Eqs. (O 
and ([4]) that additional correlation functions, i.e., GfJw), 
G^(w), and r?-(w), are introduced by the EOM method 
and therefore additional equations are required in order 
to calculate the Green's functions. First, one has to ac- 
count for correlation functions between electrons in the 
leads and the QD's, or between electrons from different 
QD's. Consider first the correlation function between 
electrons in the leads and electrons in an arbitrary QD, 



^ki( w ) =<< c ka',a} a >>■ Based on the general equation 
we find 



( W - £k )Gk»-^i^ = 



(5) 



In a similar way the general equation for the correlation 
function between electrons in an arbitrary quantum dot 
and electrons from the leads, Gf u —« a ia \c\^ a >>, can 
be obtained as 



J 



(w - Si) Gf k (u;) - - ]T V l ^Gl,^)5 il - Q(N - 1 - i)ti,H-iGSn,k(w) - Q(i - 2)t i _ 1 , i Gf_ lik (u;) = . (6) 



k' 



where r^ k (w) is another higher order correlation function. In addition we have to account for electronic correlations 
between electrons in two different arbitrary quantum dots, Gfj (<j) =<< aj CT ; aj CT >>, a correlation function which can 
be obtained from the following general equation 

(w - e 4 ) G&(u;) - U^u) - £ VkiG^a - Q(N - 1 - 2 )i M+1 Gf +lj » - 0(i - 2)t i _ M G?_ 1)i (w) = . (7) 



Note that for the particular case i — j, Eq. ([7J is identical to Eq. (0| as expected. 

Second, one has to account for additional higher order correlation functions generated by the EOM method. There 
are two different such functions, namely, r^ (w) =<< rii- a ai a ; a}- a >> and Tf k (w) =<< rjj_ CT aj CT ; c kcr >>. Applying 
the general formalism of the EOM method one can obtained two general equations for the calculation of these 
correlation functions 



(W - Ei ~ Ui) r£(w) - ^1^1 
k 

-G(N-l-i)ti, i+1 [((ni- a 
-Q(i - 2)t 4 _i, 4 



«l-crCka; a 



ai_ CT Ck-a-aio.; aj CT 



and 



(a;-e i -C/ i )rf k (cj)-^V kl ^ 1 



^'1 — crCkcrj c ]kcr 



aI_ CT Ck- CT ai CT ;c k(T \ 



c k- ff a i-^ a i^; a ]a 



(8) 



-Q(N — 1 ^/ni-o-Oi+ia;^^ + ^/at_ o .a i+1 _ .a io .; c kcr 



-0(i - 2)ti_i,i 



a i+l-cr a «-0' a jcr; C k( . 



= 0(9) 



As expected, the EOM method introduces additional higher order correlation functions. Basically, one can 
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continue to apply the same method for the calculation of 
all higher order correlation functions, with the hope that 
the resulting chain of equations will close at a certain 
level and therefore a solution for the system's Green's 
functions can be obtained. In practice it is well known 
that such a procedure will lead to an infinite number of 
self consistent equations, and to obtain a set of solutions 
for the system's Green's function appropriate approxima- 
tions has to be done. 



III. APPROXIMATION METHODS 

Several different approximations were used in connec- 
tion with the EOM method. From these, probably the 
simplest one is the Hartree-Fock approximation. If used 
for the Anderson impurity model, the Hartree-Fock ap- 
proximation simplifies the on-site Coulomb interaction 
term by replacing the number of particles operators with 
their average values, Ui < n.; CT >< rii- a >. A simi- 
lar approximation was used by Hewson to describe the 
localized magnetic states in metals (2lj. Lacroix used a 
higher order approximation to explain the Kondo effect 
in the impurity Anderson model [221 ] . 



A. Hartree-Fock approximation 

One of the most common approximation used in many- 
body theory of fermionic systems is the so called Hartree- 
Fock approximation. This approximation was used suc- 
cessfully by Hewson 2l| to investigate the single impurity 



Anderson model. Let us analyze the higher order corre- 
lation functions generated in the right-hand-side of Eqs. 
p]) and ([9]), as several of these terms are neglected in 
the Hartree-Fock approximation. The first category of 
terms we encounter are the so called normal scattering 
correlation functions. In the case of the TfJu) correla- 
tion function these terms are << 7i$— o-Qkg-; at. >> $ii 

and << ni- a ai a ; a^ a >> (l=i — 1 or l=i + 1). The 
first term describes the scattering of the conduction elec- 
trons by the localized electrons in the detector dot with 
the generation of a localized electron in the QD i = 2. 
The second term describes the scattering of a localized 
electron in the QD i — 1 (i + 1) by the electrons in 
the adjacent dot i, with a generation of an electron in 
the QD i + 1 (i — 1). In the Hartree-Fock approx- 
imation these correlation functions are replaced by a 
product of the average occupation number and an ad- 
ditional correlation function < 



and < 



ni- a ><< Ck a ;a' ja » 



n-i-a >« aia^a'ja >>, respectively. The re- 
maining terms are neglected in the Hartree-Fock approx- 
imation. A similar approximation is used for r^ k (u;). 
For this correlation function the normal scattering terms 

<< n,_ CT Ck CT ; 4*r >> and << n i-cr a lai c ka » Q =i - 1 

or l=i + 1) are replaced by < rii^ a ><< Ck^c^ >> 
and < rii-cr ><< aia'jC^ >>, respectively. Again, one 
neglects all the remaining terms in the Hartree-Fock ap- 
proximation. Accordingly, the two higher order correla- 
tion functions can be expressed in terms of lower order 
correlation functions 



$ij + EkVkiG^Sa + e(N-l- i)t i>i+1 G? +hj (u) + B(i - 2)Gf_ 



Ul — Ei- Ui 



(10) 



and 



IT k M = (n, 



Ek VkiGgkMfa + 6(JV - 1 - i)* M+ iGf +lik M + e(i - 2)Gf_ 1>k (aQ 

U3 — £i— Ui 



(11) 



We can use these approximated higher order correlation functions along with Eqs. ((3]) - (0 to calculate the Green's 
functions for the localized electrons in the system's QD's. After some simple algebra, a general recurrence relation 
can be obtained in the form 

Ai(u})G^(u) - 5 a T^(u)G^(u) - @(N - 1 - i)U, 1+1 G° +1 Juj) - Q(i - 2:./, , , = % , (12) 

I 



with 



and 



(u) - £j) {u — Ej — Ui) 
W - £i - Ui (1 - {rii-a)) 



(14) 



E 



IV, 



lk 



W - £k 



The above equation gives us a direct relation between 
fl3l the electronic Green's functions GfAcv), Gf +1 and 
Gr_r». 



5 



Let us first consider the Green's function for the de- 
tector dot, Gl^ui). First we set i = 1 and j = 1 in Eq. 
d 

- 33° ( w )] G^H - i 12 G^ = 1 . (15) 

To find the general form of the interdot Green's function 
G21 we can use the complete chain of equations generated 
for arbitrary i (i = 2, . . . , N) and j = 1. Based on this 
procedure one finds 



MP 1 - Ai(w)-Eji(o;) 



'-12 



A 2 (w) 



,(16) 



an equation which is in complete agreement with previous 
results obtained in the Hartree-Fock approximation for 
systems with one or two dots [2l|, [23| . In a similar way, 
for the final dot of the system (i — N) we find 



= A N (u)- 



t 2 

L N-1.N 



A 



N-l 



M- 



V N -2,N -1 



(17) 



Finally, let us consider the general case of an intermediate 
QD, 1 < j < N. Following a similar procedure, in the 
Hartree-Fock approximation we find 



t 



Aj +1 (u)- 



sider a set of self-consistent equations as the general fac- 
tor Ai(ui) depends on the average occupancy of the z'th 
QD, (ni_ CT ) (see Eq. (fT4|)). a value which can be evalu- 
ated using the general relation 



(18) 



A 3 M- 



Al( w )-E"(u) 



Although the general equation for the j'th QD Green's 
function looks relatively simple, in fact we have to con- 



/(w)ImG^' T (w)dw 



(19) 



f(uj) being the standard Fermi-Dirac function. More- 
over, one can see that the two spin channels, a and — a, 
are coupled in Eq. (|18[) . Such a property makes the 
self-consistent set of equations difficult to solve exactly. 
The easiest way to solve this problem is to start from 
the non-interacting case (£7, = 0), situation in which the 
set of equations fully decouples, and an exact solution is 
possible. For example, such solutions were obtained for 
the case of double or triple quantum dot systems [23|, [24| . 
The same results can be obtained using our general for- 
malism. The quantum dots Green's functions for the 
non-interacting case can be used as a starting point in 
a numerical self-consistent calculation for the interacting 
case. 



B. Beyond Hartree-Fock approximation 

The Hartree-Fock approximation lacks to explain 
many physical properties of the system. For example, in 
the case of the Anderson impurity model, Lacroix proved 
that the understanding of the_ Kondo effect requires a 
higher order approximation 22]. To obtain results be- 
yond the Hartree-Fock approximation, the terms in the 
right hand sides of Eqs. ((5J) and ([S]) have to be evaluated 
using the EOM method. The calculations are lengthly, 
but not extremely complicated, the only downside being 
that the EOM method will introduce even higher order 
correlation functions which thereafter have to be approx- 
imated. The standard approach is to neglect all higher 
correlation terms, i.e., terms which introduce correlations 
between more that four fermionic operators. 

Let us consider for example the equation for the corre- 
lation function TfJuj). Eq. © introduces an additional 
nine higher order correlation functions. We will have 
to use the EOM method to evaluate and thereafter to 
approximate these functions. The general equations for 
these correlation functions are presented in the Appendix 
IA1 Accordingly, one finds: 
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% + 5 l2 Q^ - 2)^ M Ek '^g'^ + Q(N-l- i)Ji^.S i+lij + 9(i - 2)^-^-!,,- 
'' u-ei-lh- 8 a [2X3° ( W ) + Ei(w)] - 0(iV- 1 - i)t? )<+1 M i)<+ i(w) - 0(i - 2)t?_ M M w _i(w) 
e( jV _ 2 _ z) ^±l Gf+2j(w) + e(i - 3)|E§^Gf_ 2 » 



- e< - U, - 5 a [2£°(w) + ££(w)] - 0(iV - 1 - z)if l+1 M Vi+1 ( w ) - 0(i - 2)t?_ lii M M _i(w) 
*iiE£(w) + 0(iV - 1 - i)** (n i+ i_ ff ) P M+ i(w) - Q(i - 2)t 2 i _ li (m-i-e) fy-iM 



Si-Ui- 5 a [2E£(w) + E£(w)] - 0(iV - 1 - i)t^ +1 M <ii+1 (a;) - 0(i - 2)t?_ 1 .M M _i(w) 



fa^-M + e(jv - i - z)t M+1 ^-(^) + e(i - 2)t t _ M cgH 

- C/ j; - 5a [2£°» + E£(w)] - 0(iV - 1 - z)t| l+1 M M+1 ( w ) - 0(i - 2)tf_ M M M _ 1 (u;) 



(20) 



A similar expression can be obtained also for r^ k (cj). The above equation can be combined with Eqs. ([5]) and ([7]) to 
obtain a generalized recurrence relation for the electronic inter-dot Green's function: 



LJ - £j + Ui- 



5 i1 ^1(oj) + 9(JV - 1 - Qgjjj (n i+ i- g ) (hp - 9(i - 2)^j (n»-i-<r) ViH 

- 6i - Ui - 5a [2£°(u;) + ££(w)] - 0(JV - 1 - i)i? i+1 Mi )<+ i(w) - 0(t - 2)t?_ 1>i Ar,- fi _i(u;) 



y u-ei-U- 5a [2££(w) + Si(w)] - 6(A^ - 1 - i)t? i+1 M M+1 (w) - 9(i - 2)t i 2 _ M M iji _i(w) ' ' 1 

-0(JV - 1 - i)f M+1 Gf +1) » - 9(i - 2)i,_ M G?_i,» 



G?» 



0(jV_2-i)|±i^^<n 4 - cr ) 



E ( - Ui - 5a [2£°» + E£(w)] - 0(JV - 1 - i)t? i+1 M M+1 (a;) - 0(< - 2)t?_ 1 ,%iH 



w - e< - Ui - «Sii [2££(w) + E£(w)] - 0(JV - 1 - i)^ ii+1 M ili+ i(w) - 0(< - 2)t?_ M M<,j_i(w) 
fa^M + 0(JV - 1 - O^iB^fw) + 0(i - 2)t,_i,,C§(a;) 



-u 



w - 6j - ^ - 5 a [2£°(w) + Ei(w)] - 0(iV - 1 - i)tf )i+1 M i!i+1 (a;) - 0(< - 2)t?_ M M M _i(w) 

Ui (ni- ff ) 



" ' ' w - £i - Ui - <5a [2£° (w) + £* (w)] - 0(AT - 1 - ^ 



- £i - Ui - (5a [2£°H + Ei - 0(7V - 1 - £)t? i+1 M M+1 (w) - 0(i - ^^AnH 



(21) 



The above equation in its most general form it is a self- 
consistent integral equation. Different from the Hartree- 
Fock approximation case, here, the inter-dot electronic 
Green's function GL(uj) depends on the inter-dot elec- 
tronic Green's functions for the near neighbors G? ±1 Jlj) 
and the inter-dot electronic Green's functions for the 
next-near neighbors Gf ±2 j (w). Similar to the Hartree- 
Fock approximation, the inter-dot electronic Green's 
function depends on the average electronic occupancy 
(nj_ CT ) in dot i, which can be calculated using the general 
relation (|19l) . Eq. (|2ip can be used to obtain the intra- 
dot electronic Green's function Gfj(w) by simply setting 
j = i. Additionally, anomalous terms (Afj(u>), BL(w), 
and GfAu))) are present in the general equation for the 
electronic Green's function Eq. (|2Tj) . These terms will 

introduce new average values of the type (aj_ CT Ck_ 



c^._ tr (ii—(T /) an d (a|_„a,'- CT ), which have to be calcu- 



lated self-consistently. For example, one can calculate 



a\_ a c k - CT 



or ( a- „ai-,T ) as 



/(w)ImGfc> 



(22) 



/(w)ImG- ff (w)da; . (23) 



Although neglected in some quantum dot systems studies 
[Iol | , such anomalous averages have proven to be of main 
importance in the investigation of the system's physical 
properties [22|, [25j. 
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IV. T-SHAPED DOUBLE QUANTUM DOT 
SYSTEM IN THE COULOMB BLOCKADE 
REGIME 

The simplest T-shape quantum dot system described 
by the model is the double quantum dot system. Sim- 
ply, if we set t 2 3 — 0, the only remaining components 
of our system are the quantum dot connected to the ex- 
ternal leads (detector dot) and a second quantum dot 
connected to the first one, but not to the external leads 
(side dot). In the following we will investigate the elec- 
tronic density of states in the detector quantum dot 
(Pl(u>) = —hxxGxx(u})/Tr) and its dependence on the on- 
site Coulomb interaction terms U\ and U2 using various 
approximations. First, we will point to various levels 
of approximation for this system and thereafter we will 
present numerical results for the occupation number of 
the system. 

The simplest possible situation is for U\ = U2 = 0, 
case in which the density of states for the detector dot 
can be calculated analytically as: 

= ~ ^—2 , (24) 



where A is the imaginary part of + it]) (77 — > 0),i.e., 
A = 7r^|F kl | 2 <5(w-£ k ) , 

k 

and d(x) is the delta Dirac function. In this situation, 
the electronic density of states will present a double peak 
structure, the structure and the position of these peaks 
being controlled by the system's parameters (23|. Let us 
consider the case of finite on-site Coulomb interaction in 
both the detector (U\ 7^ 0) and side (U2 7^ 0) dots. 

A. Hartree-Fock Approximation 

The general theory presented in the previous section al- 
low us to estimate the detector's Green's function within 
the Hartree-Fock approximation: 



1 

(k>-ei)(u-ei-t/i) _ ,2 u-e 1 -U 2 {l-{n2-< r )) v (. ,\ 
u>-ei-(7i(l-<7U-„» 12 (w-e 2 )(w-e 3 -l/2) ^ > 

(25) 

Although the equation for the detector's dot Green's 
function is analytic, before one can estimate the detector 
dot density of states pi (uj), one has to self-consistently 
evaluate the average occupation number for both the de- 
tector and side dots. In the case of the side dot, the 
required Green's function is given as 

G^H = 

1 

(ui-e^ju-ez-Uz) tja 

(26) 

Both occupation numbers can be calculated using Eq. 
(|19p by means of an iterative numerical evaluation. The 
resulting density of states presents additional peaks as a 
result of non-zero on-site Coulomb interaction terms. 



B. Beyond Hartree-Fock Approximation 

In general, the Hartree-Fock approximation can give 
some insight on the importance of the on-site Coulomb 
interaction terms U\ and U2, however, it is necessary to 
go beyond this approximation to recover important phys- 
ical properties of the system. For example, Lacroix [22[ 
proved that a calculation of the Kondo temperature us- 
ing the single impurity Anderson model in the U\ — > 00 
limit requires terms beyond the Hartree-Fock approxima- 
tion. We expect that such terms will be relevant also to 
the finite U\ and U2 case. To calculate the detector and 
side dots Green's functions we will need to use a set of 
coupled equations: 



u - £1 - E£(w) + f/i 
ii 2 G£iM - Ui 

: 1 



££M+t?2 (na-a)i\ 2 (w) 



ui - ei - Ui - [2££H + E£(w)] - t 2 12 M 12 (ui) 
w - £1 - U ± - [2E°(cj) + E{ {u})} - t 2 12 M 12 (oj) 



Ui (rii_ CT ) 



u-ei-Ut- [2X3^(w) + Ei(w)] - t 2 12 M 12 (uj) 



(27) 
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and 

V w - e 2 - (7 2 - q 2 M 2i (w / 

| fig g^Nga (ng-g) _ fia^H 
w - ei w - e 2 - U" 2 - if 2 M 2 i(w) uj - e 2 - U 2 - t\ 2 M 21 (uj) 

= 1 , U2 {n 2 - a ) 

u — e 2 - U 2 - tj 2 M 21 (to) 



Clearly, the above equations are self-consistent due to 
the presence of terms involving the occupation numbers 
in the detector and side dots. Moreover, the presence 
of additional correlation functions, Gf 2 (w) and G 21 (w), 
makes the calculation even more complicated. Additional 
equations are required to replace GJ 2 (w) and G 21 (w) in 
terms of G^uj) and G 22 (u;). After some algebraic ma- 
nipulations, both detector and side dots Green's func- 
tions can be expressed only in terms of the average occu- 
pation numbers (ni CT ) and (n 2<7 ) and a series of anoma- 
lous averages of the type (a\_ rT c\ f .- (T V (^i_ a Ck-a), 

(^a2_ a Ck-o-\ and (a^^ai-cr^, along with their hermi- 

tian conjugates. The presence of these terms makes the 
problem self-consistent as such averages have to be calcu- 
lated using their corresponding Green's functions using 
Eqs. (HH, (EH), and (USD. 

Two levels of approximation can be considered in this 
limit. First, one can consider all the anomalous averages 
to be zero. A similar procedure was used by Zimbovskaya 
10] in connection with the Coulomb blockade regime in 
single quantum dot systems. On the other hand, Lacroix 
proved that the Kondo regime (U\ — > 00) of the sin- 
gle impurity Anderson model can be explain only if the 
anomalous averages are taken into account [22j ■ Our cal- 
culation can be reduced to the single impurity Anderson 
model if we consider the limit t\ 2 —¥ 0, i.e., we decou- 
ple the side dot of the T-shape system. The obtained 
Green's function for the detector dot matches the previ- 
ous findings of both Zimbovskaya and Lacroix, assuming 
the right levels of approximation. 

In the case of the double dot T-shape system, for fi- 
nite onsite Coulomb interaction in both the detector and 
side dots, one has to use numerical calculations to eval- 
uate the system's properties. We choose to calculate the 
system's average occupation number, n t . a ={ni a )-\-(n 2a ) . 
Our main goal is to prove the differences between the 
two levels of approximation respect to the anomalous av- 
erages mentioned above. First, we followed Zimbovskaya 
and we disregarded these averages (BHF1). Second, we 
followed Lacroix and we took into account these averages 
(BHF2), however, in this limit we considered a more gen- 
eral situation with finite on-site Coulomb interaction in 
both component dots (not only the particular situation 
with infinite on-site Coulomb intercation). 



(28) 

I 

Figure [5] showcases a numerical estimation of the elec- 
tron occupation number per spin orientation of the dou- 
ble dot T-shape system. We considered two different sit- 
uations, A/t = 0.5 and A/t = 3, corresponding to the 
case of a slow (weak coupling) and fast (strong coupling) 
detector configurations. The plot presents the Hartree- 
Fock approximation results (dashed line) and the results 
obtained using terms beyond the Hartree-Fock approxi- 
mation without - similar to the approximation in Ref. [l^ 
- (dotted line) and with (dot-dashed line) considering all 
anomalous averages in the detector and side dots Green's 
function. For comparison, we plotted also the simple 
case without on-site Coulomb interaction, U\ = U 2 =0 
(dot-dashed line) , situation in which the Coulomb block- 
ade regime is not possible. Clearly, the inclusion of the 
anomalous terms in the calculation resulted in significant 
changes in the system's total average occupation num- 
ber, especially in the Coulomb blockade regime. When 




-20 -15 -10 -5 5 10 -20 -15 -10 -5 51 



FIG. 2: The system's average total occupation number (per 
spin orientation) as function of the localized energy level in 
the detector dot in various approximations: Hartree-Fock ap- 
proximation (dashed line), beyond Hartree-Fock approxima- 
tion without (dotted line) and with (dot-dashed line) consid- 
ering all anomalous averages. The considered parameters are 
Ui/t = 10, U 2 /t = 5, E F /t = 0, T/t = 1(T 4 , A/t = 0.5 
(left), and A/t = 3.0 (right). We also considered Ei/t = e%/t. 
For comparison, we present results in the simple case without 
on-site Coulomb interaction Ui = U2 = (dot-dashed line). 
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the energy of the localized level in the detector dot is 
well bellow the Fermi energy of the free electrons on the 
leads (Ep = 0), one electron per spin orientation will be 
accommodated on each localized level of the system. As 
the energy increases, the effects of the on-site Coulomb 
interaction are becoming important, so only one electron 
can occupy the localized levels. Once the energy of the 
localized level raises above the Fermi energy, the occupa- 
tion of the localized levels approaches zero as electrons 
will flow directly from one lead to the other. In princi- 
ple, the relative value of the localized level characteristic 
energy can be controlled using an external gate. 

V. CONCLUSIONS 

In conclusion, we presented a detailed analysis of a T- 
shaped multi-quantum dot systems using the equation 
of motion method. All our calculations are done in the 
case of a finite on-site Coulomb interaction in each of the 
component dots. The results are obtained within various 
approximations and for each case we provided a set of 
self-consistent equations which allow the calculation of 
the electronic Green's function in the system's compo- 
nent dots. The simplest approximation is similar to the 
one introduced for a single quantum dot system by Hew- 
son [2l| and is equivalent to the Hartree-Fock approxima- 
tion. In this case, the general recurrence equation which 
relates the electronic Green's function between adjacent 
quantum dots in the system can be solved and we ob- 
tained a set of self-consistent equations for the electronic 
Green's functions. The self-consistency in the case of 
the Hartree-Fock approximation is introduced through- 
out the electronic occupation number in each of the com- 



ponent dots. When an approximation which includes 
terms beyond the Hartree-Fock approximation is consid- 
ered, the situation is more complex and the analysis of 
the results is more complicated. Even in this case, our 
calculation lead to a recurrence relation between the elec- 
tronic Green's functions of the system's quantum dots, 
however, different than in the Hartree-Fock approxima- 
tion, the electronic Green's function for a component dot 
relates to the electronic Green's functions in the near and 
next-near quantum dots. The main result of our calcula- 
tion is the inclusion of several anomalous average terms, 
terms without which important physical effects such as 
the Kondo effect cannot be evaluated properly. 

As an example for our general theory, we considered 
the case of a double-dot T-shape system consisting on a 
detector dot connected directly to the external leads and 
a side dot connected only to the detector dot and not to 
the external leads. In both system's dots we considered 
a finite on-site Coulomb interaction and we discussed in 
various approximations the total average electronic oc- 
cupation of the system's localized levels. We found that 
the inclusion of the anomalous terms in the evaluation of 
the detector and side dots Green's functions is very im- 
portant in the evaluation of the system occupation num- 
ber. Depending on the relative position of the detector 
dot localized level respect to the Fermi energy, different 
occupation regimes are possible, pointing the important 
role of the on-site Coulomb interaction term. 
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Appendix A: Higher order correlation functions 



In this Appendix we present the general procedure for the estimation of the four particle correlation function 
Tfj((j) =« ni- a a,ia\cvA a » (see Eq. (JSj) ) using the EOM method beyond the Hartree-Fock approximation. A 
similar calculation can be performed for the correlation functions in Eq. The general equations for these 

functions are 



(w - e k ) ((ni-o-Cko-jat^)) = Vki ((ni_ ff ai CT ;a 



+ Y>k'i 



k' 



ai-o-Cko-; a\ a 



+t 



12 



[((»! 



(Al) 



(cj - £k) ((aJ_ .Ck-ffa lo -;ajj (r 



al-crCk-o-) Sij + Vki ((ni-o-ai CT ; a] 



k' 



V V1 

[«"! 



r, 



(A2) 



and 



(to + s k - 2ei - Ui) ((cl._ a a l - (T ai rT ;a 1 



3<? 



"k-cr a i-W hi - Vki ((ni-o-aiaja 1. 



ci-a C V-<r a lv'> a ]<j 



+h2 [((4_ CT ai_ CT a 2CT ; 



For 1< i < iV - 1 



+U,i+1 

+Q(i - 2)U- lti 



a\_ a a i+ i- a a i+ i a ] a^)) + ((ni- a a ia ;a^ ja . 



a !+l-o- a i-o- a i+lo-; a] CT 



al_ <T a i - 1 - a a i+la ;a] ja 



a I-l-o- a i-o- fl i+lo-; a j<j 



+9(iV - 2 - i)ti+i,i+2 ( (ni- a a i+ 2a\ a 



3 a 



(u - e i+ i + Ui) 



al_ CT a 2 -aCk' CT ;a} \\ - ((4,_ iT a2- (7 ai CT ; aj ff 



+e(i-2)tj_i, 

+6(iV - 2 - i)i i+M+2 ((at_ CT ai + 2- CT aiCT; at 



ai-<T a »+i-<T a »-i^; a J<7 



and 



(cj - 2e 4 + £j + i - Uj) 



4+l-o- a i-o- a io-i a lcr 



a l+l-cr a i-cr/ % 



a i+l- ( r a i-cr a i+lo'; a]cr 



+e(z-2)t l _ 1 

-6(7V - 2 - i)t i+1)i+2 ((aj +2 _ CT a i _ (T a i(T ;a] 



a i+l-cr a i-o- C k:'<T; a] CT 



a i+l-a C k' -crO-itr', 



For 2 < i < iV we have 

(w-e 4 _i) ^n J _ CT a J _i cr ;a] c 



+U-i 

+0(N - 1 - i)U,n-i 

- 3)ij_2,,-i Uni-aai^; a j CT )) ' 



a l-i-<T a i-<T a i-i<x; a ] CT )/ + ( (ai- ff a i-i-T a *-i^; a ] 
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J" 



a i-a a i-l-<r a i-l-<?> a ]a)) + ( \ni-<T a i(r\ a j 



Tli—l — cr&icri ® i 



and 



+e(jv-i-»)t i)i+ i 

+0(i - 3)ti_ 2 ,i_i ^/aj_ CT ai 



a i-<T a i-l-o- a i+lo-; Ojo- 



a i+ i_ <T a i _i_ cr ( 



-2-a-CLicr] 0. 



(A8) 











Oi-l-o 









Si2 ^2 Vu>i \( c k'-o- a ' ; -^ a ^; a 



a-i-i-o-ai-crai-icr] a 



i)) + ((^-i- CT a«x;<4)) 



+e(jv - 1 - »)t 4fi+ i 

-6(i - 3)ti_ 2 ,i-i ((aj_ 2 _ .ai_ -a i0 .;aj 



ai-l-o-Oi+l-o-aicr; 1 



(A9) 



In all the above equations higher order correlation functions involving six particle processes were neglected. Obviously, 
to truncate the chain of general equations produced by the EOM method one has to make additional approximations 
in Eqs. (|Alj) - (|A9j) . In the following, the standard procedure we used to approximate the four operators correlation 
functions, << A^a-B- a C a ; D a >>, is to replace them by a product of the average value of two operators and the 
correlation function of the remaining two, < A_ (T _B_ . ><< C a ;D a >>. Of course, this procedure is not applied 
when any of the four particle correlation functions resembles one of the initial functions we try to estimate, for the 
particular case of Eq. (|8]), TfJtu) =<< rt^o-a^; a. >>. Based on this approximation we rewrite Eq. 



as 



[u-ei-Ui- S tl (2E£(w) + Ei(w)) - 0(JV- 1 - i)<«%i(w) - 9(* - 2)i?_ lii M < , i _ 1 (w)] r?.(w) 



(ni-a) 



Sij + 5 i2 Q(i - 2)ti_i,. 



ui - e. 



+ Q(N-l-i) 



i.i+l 



+ 0(i - 2) 



CJ - £i_i 



+(-)(. V -2-i)-^±^Gf +2> ,H + e(z-3)^^G; , 



w - e-i+i " w — £j_i 

-G£(w) [fa5£(w) +6(iV - 1 - i)t? (<+1 (m+i-a) Pi )i+ i(w) - 0(i - 2)*?_ M (ni_i_ ff ) P<,i-i(w)] 
+ 9(iV - 1 - i^i+xBf^u;) + Q(i - 2)ti_i )< C^(w) , 

where for simplicity we introduced the notations 

My(w) 



1 



1 



1 



and 



CJ — £j CJ — £j + Ui UJ — 2Ei + Ej — Ui 

1 1 



u — Ej + Ui uj — 2e-i + Ej — Ui 



(A10) 
(All) 

(A12) 



k 



w - £k 



cj + e k - 2ei - C/i 



a 2-cr c k-(T 



uj + £k - 2ei - C/i 



Ek' ^k'l (<4_ CT Ck'- CT ) - 



c{.,_ a a X - a 



UJ - £ k 



w - £k 



a 2- CT a l-cr , 



(A13) 
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and 



B^M = 8 rj + 5 n J2 VviGLtb*) + 0(i - 2)t i _ 1)i Gf_ 2J (a;) 



a j+l-cr a *-o- , 



w - Ej+i + Ui lo -2si + e l+1 - U, 



'k'-CT u *-o- 



+ 0(i - 2)t 4 _i,i 



'i-l-cr^-o- , 



(^i-er^+l-er) f 1 

\ / Vw - e i+ i w - e l+ 

— U,i+1 \ a l+l-a^- 



- e i+ i + J/, 
1 



4'-<7°2- 



w - e l+ i ui - 2ei + e l+ i - Ui 



a 2-cr C k'-(T 



ui - e i+ i + t/j w - 2ei + e l+1 - U, 



a i+l-cr a i-l-0' 



w - e 4+ i + C/i a; — 2e. t + e i+ i - U t 



+ t i+ i, i+2 e(N - 2 - i) 



a i _ .a l +2- ( r 



l i+2-(T u '*-o- 



w - £i + i + Z7j a; - 2e 4 + e l+1 - Ui 



ck(w) = (% + e(jv - 1 - t)*i,i+i) 



'i-l-cr^-^ 



UJ - £j_i + J7j CJ - 2£j + £i_i - C/j 



4_ CT c k '_ c 



'k'-o 



k'l 



cj - £j_x + [7j w - 2e l + £j_x - Z7i 



-6(AT-l-i)t 



a 7:+i-r7 a i-i- 



a i-l-<r a i+l- 



»,»+! 



W - £j_l + C/i W - 2e,; + £ i+1 - fjj 



+e(t - 3)ti_ 2 ,i_i 



4-<r a i-2- 



a i-2~a a "t-o 



ti—i a ( a.- „a,:_i 



CJ-£i_l + Uj w - 2e, + £,_1 + Ui 
1 1 



■ il - 1 ''( a «- 1 - fl -')(^ + 



UJ - £,_i W - £,_1 + 
1 



+Q(N- !-«)■ 



W - 2E; + £,_! - [/j 
t 



W - £i-l 



1-CT 



i+l-cr a »-o- 



(A14) 



(A15) 



Additionally, similar to Lacroix 22j we introduced two terms related to the interaction between free electrons in 
the leads and the localized electrons in the detector quantum dot: 



UJ + £k - 2£i - Ui 



(A16) 



kl 



Ek' ( C L- ff C k-<r) Ek' ( C L<x C k' 



W - E k 



(A17) 



As a general remark all the above terms introduce anomalous correlations between electrons in different quantum dots 
or between electrons in the leads and in the detector quantum dot. Such terms are usually neglected; however, it was 
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proved that they are of main interest when physics associated with the Kondo effect is discussed [22]. Consider for 
example the simple situation in which the T-shape system consists only in a single quantum dot (ti2 — 0) and focus 

on the anomalous terms introduced by Eq. (|A17[) . i.e., (c£/_ CT Ck-<7^ and /cj c _ .ci c /_ -^. The simplest approximation 

will be to replace these terms with <5kk'/(£k), although a more careful analysis shows that 



with 



Ck- 



Jkk' 



f(u) Im ( ( c k - a ; cj. 



VklHc' 



w-e k (w - £k)(w - £ k ') 
Under this approximation, if we consider terms up to the quadratic order in Vki we find 



A 

n 



1 



w - 2ei - t/i 
2vriT 









HG- 




— 27T 




2mTj 





(A18) 



(A19) 



(A20) 
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